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An extension of the notion of concurrence introduced by Wooters is used to quantify the entangle- 
ment of certain multipartite pure states, namely, the BCS state of superconducting compounds. This 
leads to a definition of the macrocanonical entanglement of pairing (MEP) for which we compute 
an analytical formula in terms of two adimensional numbers, the cut-off and gap numbers, which 
depend on measurable physical quantities. We find that strongly coupled BCS elements like Pb and 
Nb have much larger MEP values than more conventional BCS transition metal superconductors. 
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The study of entanglement of quantum states is now a 
central issue in the modern theory of quantum informa- 
tion Entanglement is not seen just as a sort of pecu- 
liar or curious property of multipartite quantum states, 
but more importantly, it has become a resource to achieve 
novel tasks, such as teleportation, quantum crytogra- 
phy and other quantum transmission protocols surpass- 
ing the capabilities exhibited by their classical counter- 
parts. Moreover, entanglement is also something of direct 
experimental importance since it amounts to interaction 
among parties: just to mention a simple instance of this, 
the application of the CNOT gate to a factorizable state 
■= 2-i/2[|o)_|i)]|i) produces oneofthe Bell states, 
namely, t/cNOT|*)i„ = 2-i/2[|01) - |10)] := I*"). In this 
process, the output state is more difficult (and expen- 
sive) to realize than the initial state |^)in because of the 
CNOT gate that implements entanglement. Therefore, 
to know wheather a given purported state is entangled 
or not, as well as how much entangled it is, is of great 
importance both theoretically and experimentally. 

In condensed matther theory, we are used to deal 
with quantum many-body states (or multiqubit states) 
in which their strong quantum correlations are responsi- 
ble for novel properties or states of matter, like quantum 
liquids, or quantum phase transitions Example of 

strongly correlated states abound in these areas, for ex- 
ample, valence bond states are nothing but Bell states, 
and properties like entanglement swapping correspond to 
resonating valence bond (RVB) state configurations. 

In low dimensional systems, like quantum spin chains 
and ladders, it is known that the effect of quantum fluctu- 
ations is stronger than in higher dimensions: the factor- 
izable Neel state is a good starting point to describe the 
ground state of the antiferromagnetic Hcisenberg model 
in D = 2 or more dimensions, but it is unsuitable for 
D = 1 where the Bethe ansatz solution is a complicated 
superposition (entangled) of single particle states. Thus, 
it is natural to ask whether the new ideas about qual- 
ifying and quantifying entanglement that have emerged 
in the field of quantum information can be helpful to de- 
scribe the complicated patterns of behaviour exhibited 
by strongly correlated systems in condensed matter. 



Recently, the entanglement properties of the one di- 
mensional XY model in a transverse magnetic field have 
been analyzed in the vicinity of a quantum phase transi- 
tion by Osterloh et al. ^ and Osborne and Nielsen [||. 
In this model, entanglement shows scaling behaviour near 
the transition point and remains short ranged 0. The 
quantification of entanglement is made with the entangle- 
ment measure known as entanglement of formation E{p), 
introduced by Bennett et al. § to describe the resources 
needed to create a given entangled bipartite state, either 
pure or mixed p [||. 

Generally, it is difficult to find closed mathematical 
expressions of E{p) solely in terms of p, but for the special 
case of mixed states of bipartite qubit systems Wooters 
§, § found one such a formula. This formula makes 
use of what Wooters calls a spin flip transformation, 
defined as 

\^):=ay\r) (1) 

where cr., := f „ | is the Pauli matrix in the com- 

putational basis {|0) := It)^!) '■= li)}- For a gen- 
eral state p of two qubits, the spin-flipped state is 
p ■= {(Jy (8) ay)p*{ay ay). When the state is pure 
p = the entanglement of formation can be written 

as E{p) = £{C{ip)), where the concurrence C is defined 
as 

C{4'):=\{m\ (2) 

&nd£{C) := /i(i[l-hVl-C2]), h{x) := -xlog2a;-(l- 
x) log2(l — x), is a monotonically increasing function of 
C that ranges from to 1 as the concurrence goes from 
to 1. 

In fact, Wooters proposes to use concurrence (||) as an 
entanglement measure in its own right. Here, we shall 
adhere to this proposal by studying the concurrence of 
a physically realizable state such as the BCS state [ pl)[ . 
There are very few cases where we have a solution to a 
quantum many-body problem in the form of an explicit 
wave function. The BCS theory of standard supercon- 
ductors provides us with one of these examples. Specif- 
ically, the general quantum state representing a super- 
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conductor carrying a supercurrent at T=0 temperature 
is 

k 

where |0) denotes the zero-particle Fock state. In this 
state the electrons are created in Cooper pairs for the oc- 
cupied states with quantum numbers (k+Q, "f; — k+Q, I 
). All these pairs have the same momentum 2hQ. This 
pair momentum represents the finite supercurrent and it 
is usually very small. We shall concentrate in the static 
condensate of Cooper pairs with zero supercurrent Q = 0. 
The parameters Uk and Vk represent the probability am- 
plitudes of creating quasi-holcs and quasi-electrons, re- 
spectively. They satisfy the following properties: 

Uk,vk e M; fc |k| 

The first condition comes from the normalization of the 
state |BCS)6(, and they only depend on the modulus of 
k. The phase factor c'^ is arbitrary, but is the same 
for all Cooper pairs. In the macrocanonical BCS state 
the number of Cooper pairs N is not a well defined 
quantity. By series expansion, the state can be thought 
of as an average over an ensemble of states |iV, Q) with a 
definite number TV and pair momentum 2hQ: |BCS)0 := 
j:j,e'^<^AN\N,Q),j:j,A% = l. 

The advantage of having the explicit form of the 
ground state wave function is that we can compute any 
quantity needed for an entanglement measure. In partic- 
ular, it is possible to compute the reduced density matrix 
/9(ki,kj) by tracing out over all Cooper pairs with mo- 
menta k ^ ki, kj. This is a bipartite density matrix for 
which explicit formulas for the entanglement of formation 
also exist in terms of the concurrence 1^. However, we 
notice that in the case of the macrocanonical BCS state, 
the matrix p(ki,kj) corresponds precisely to the state 
formed by the product of two Cooper pairs with momenta 
(ki,kj), namely, [ufejO) -t-Ufc; «) [w/cjO) -|-UfeJl)], where 
here |0), |1) denote states with zero and one Cooper pair, 
respectively. Then, the reduced density matrix leads to 
the same original problem we are dealing with, but with 
only two pairs. 

Despite of this difficulty, is it still possible to use the 
concurrence to devise an entanglement probe for the 
macrocanonical BCS state? In principle it looks difficult 
since the BCS state is a many-body (multiqubit) state 
and it is known that concurrence fails to capture entan- 
glement properties of multiqubit states. For example, it 
is known that any qutrit state can be entangled in two dif- 
ferent ways |l| : cither as a I GHZ) := 2-i/2(|000) + |lll)) 
state or as a Werner |W) := 3-^^^{\001) + |010) + |100)) 
state. Both of them yield zero concurrence since |GHZ) = 
12-1/2(1000) - |111)) and |W) = i3-i/2(|iio) + |101) + 



|011)). Thus, concurrence docs not detect the existence 
of entanglement for qutrits. 

However, here we show that relying on physical 
grounds and motivated by the physical meaning of the 
concurrence, it is possible to give an entanglement mea- 
sure for the BCS state based on the notion of concur- 
rence. We do this in two steps. 

In what follows, we shall study a many-body state |4') 
which is a BCS ground state with 9^0: 

Ivf) := |BCS)o = HK + yAA^Ol^) (5) 

k 

This BCS state is the solution of minimum energy to a 
reduced Hamiltonian called pairing Hamiltonian [|2| : 

ifrcd J2 ^^kblbk - J2 ^k.k'&L^k (6) 

k k^k' 

where b^^ I'^-k j.' '■= c-k,iCk,| are operators cre- 

ating and annhilating Cooper pairs, respectively. The 
solution to this variational problem yields the following 
expressions for the probability amplitudes 

(7) 

Ek = ^Jel + Al 

where Ek is the energy of the quasi-particles (excita- 
tions), and Afe is called the gap function, which is deter- 
mined by the self-consistent solution of the gap equation 
Afc = ^ X^k' ^T^.k'- This solution represents a BCS 
superconductor or SC state. 

In the first step, we extend the notion of concurrence 
to many-body states based on the physical interpretation 
of concurrence introduced by Woofers . Namely, for a 
spin i particle the spin-flip operation (|l|) is the time- 
reversal operation. For two-qubit states, we can argue 
that concurrence can serve as an entanglement measure- 
ment directly form this definition, without resorting to 
its connection to £{C). The rationale goes as follows: as 
lip) is obtained from ji/)) by time inversion, we intuitively 
expect that if {ip) is very much entangled, then lip) will 
be very similar to |-0) thereby C{ip) ~ 1. On the con- 
trary, if I?/)) is factorized into two states then 1-0) will be 
very different from \ip) and then C{ip) ~ 0. When the 
given state is very entangled, the timc-rcvesed state is 
very close to the original state and their overlap is very 
large. 

Thus, it is reasonable to extend the notion of spin- flip 
opcrtation by the time-reversal operation. Let us define 
the concurrence C(BCS) of the BCS state (|) by means 
of the overlapping with its time-reversed state, namely, 

C(BCS) := |o(BCS|BCS)o| (8) 
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with 




1 


|BCS)o Ut\BCS) 


(9) 
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The action of the time-reversal operator Ut on position, 




momentum and spin variables is |]l3| 




BC 


UrriUl = Tj 






UTk,ul = k. 


(10) 
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In the second step, we realize that in the absence of 
a precise connection between entanglement of formation 
and concurrence beyond two-qubit states as the one pro- 
vided by Woofers §, 0, wc must define an entangle- 
ment measure with respect to a reference state for which 
we know that it has zero entanglement. The natural can- 
didate for this is the Fermi Sea state |FS) defined as 



IPS) n iriiio) 



(11) 



k<kF 



Using ([lO|), this state has maximum concurrence 
C(FS) = 1 despite being unentangled. Thus, we choose 
as our definition of entanglement what we shall call MEP 
macrocanonical entanglement of pairing E(BCS) defined 
as 



i(BCS) := (FS|FS) - (BCS|BCS) 



(12) 



With this definition the Fermi Sea state has zero entan- 
glement. Now, we compute the MEP to check that it 
serves as a record signalling and quantifying the onset 
of entanglement in the macrocanonical BCS state. Us- 
ing the fact that Ut\0) = |0) and (|o|) along with the 
canonical anticommutation relations, we find 



iBcs)o = n( 



Uk 



«fc4.TC-ka)|0) = |BCS). 



(13) 



This means that under time-reversal, the state with 9 — 
is transformed into another one with 9 = n. Let us point 
out that all states |BCS)6/ have the same energy with 
respect to the reduced Hamiltonian (||) p4[ . 

Now, we can compute the concurrence of the BCS state 
as the overlapping between states (llSh and (@) , yielding 



(14) 



C{BCS)^Y[\iul-vl)\ 



This suggests to introduce the notion of a partial con- 
currence Cfc(BCS) for each Cooper pair of momentum fc, 
and using the BCS solution (0) to the gap equation we 
readly find 



Cfc(BCS) ■.= \{ul~vl)\ 



(15) 




FIG. 1: Partial concurrences Cfc(BCS) in terms of energy 
differences w.r.t. the chemical potential /i, for several values 
of the superconducting gap A. 



In Fig. we plot Cfc(BCS) as a function of the energy 
for several values of an homogeneous gap function := 
A. It clearly shows that the partial concurrences are 
different from 1 in the vicinity of the Fermi surface, and 
this region extends within an interval of the order of 2 A. 
This is precisely the region where Cooper pairs are being 
formed. Thus, deviations of Cfe(BCS) from 1 allows us to 
detect the onset of correlations between pairs of particles. 
In fact, when the gap vanishes the solution to Uk,Vk in 
(0) yields Cfc(BCS) = l,V/c, which agrees with the fact 
that it represents a normal metal with an uncorrelated 
(factorizable) ground state (|ll|). 

We can proceed even further and compute the total 
concurrence of the BCS state. It can be computed ex- 
actly in the continuum limit (^) J d^k with the 
following result 



C(BCS) 
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ni :=N{eF)huJu 
n2 :=7V(eF)A 



(16) 



where Ni^e-p) is the density of states at the Fermi level and 
Lo-Q is the Debye frequency. For instance, for a parabolic 
dispersion relation like ek we have 



_L3 

27r^ 



, dk 
dek 



U 



mk-p 
2^ 



(17) 



Therefore, we find that the BCS concurrence (|T^) de- 
pends on two adimensional quantities with a physical 
origin that we call the cut-off number ni and the gap 
number n2. From (|l^) we see that this concurrence is 
always < 1, and the maximum value is attained when 
712 = corresponding to the absence of superconductiv- 
ity (and thus, no correlated pairs). In Fig. || wc plot 
the macrocanonical entanglement of pairing MEP (|l|) 
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to show how it depends on the numbers ni and n2. We 
see that for a fixed value of ni (and thus of the cut-off 
lud and N(eF)), the entanglement increases with 712 and 
in turn with the superconducting gap A at T = 0. 



E(BCS) 




FIG. 2: Macrocanonical entanglement of pairing (MEP) 
E(BCS) as a function of the cut-off number ni and gap num- 
ber 712 (p^. The blue regions represent MEP high while red 
regions correspond to MEP low. 



The nice feature about the analytical formula ( J16| ) for 
the entanglement is that it depends through m and 712 on 
the experimentally accessible quantities iV(eF), wd and 
A. We propose to use MEP in superconducting com- 
pounds as an indicator of the entanglement or quantum 
correlations due to the presence of Cooper pairs, ft is 
interesting to start with the simplest cases, namely, the 
superconductive elements. In Table || we have computed 
the values of MEP from available experimental data , 
| p6[ , for transition superconducting metals and lead. 
We find relevant to put the MEP values together with 
the dimensionless electron-phonon coupling constant A. 
In the conventional BCS theory, we have A ^ I and 
the effect of phonons is to provide a cut-off Hlub to the 
possible electron energies. However, there are materials 
like Pb and Nb for which retardation effects of phonons 
are relevant for the pairing of electrons, and they need 
to be treated dynamically within the Eliashberg exten- 
sion of the BCS theory |T^. In this theory A is defined as 
A := 2 doja'^ {uj)Nph{uj) /ui where iVph(ti-')is the phonon 
density of states and a^{uj) is the electron-phonon cou- 
pling strength. In the weak coupling limit A ^ f, it 
reduces to the BCS coupling parameter A « iV(eF)Vb. 

From this table a clear picture emerges: as A gets big- 
ger, there is a large increase in the value of MEP, specially 
for the case of Eliashberg superconductors which have 
values of MEP about 3 orders of magnitude higher than 
in more conventional superconductors. For the group 
of elements Ru, Mo and Os that have similar values of 
A, their correspoding MEP values are also very close. 
Thus, we conclude that strongly coupled BCS supercon- 



Superconductive Elements 


-logio(E(BCS)) 


A 


Tc(K) 


Transition SC Metals 


Hf 


7.557 


0.14 


0.13 


Ru 


6.813 


0.38 


0.49 


Mo 


6.381 


0.41 


0.92 


Os 


6.290 


0.44 


0.66 


Eliashberg SC 


Nb 


3.549 


0.82 


9.25 


Pb 


3.295 


1.55 


7.20 



TABLE I: A list of superconductive elements with their values 
of macrocanonical entanglement of pairing E(BCS), electron- 
phonon coupling constant A and critical temperature Tc- 



ductors arc characterized by large MEP values. We find 
this reasonable since the effect of phonons is to enhance 
the electronic correlations ||l^ . This enhancement is also 
responsible of the higher Tc values in Table |. 

It might be possible that MEP could also be computed 
for other more complicated superconducting compounds 
such as heavy fermion materials, high-Tc cuprates, MgB2, 
fullercnes etc., with the purpose of having an indicator 
to distinguish them in different categories. Specially in- 
teresting could be the dependence of MEP with T as 
we approach the transition temperature Tc from the SC 
phase. 
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